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Abstract 

We say that a quantum spin system is dynamically localized if the time-evolution of local 
observables satisfies a zero-velocity Lieb-Robinson bound. In terms of this definition we have 
the following main results: First, for general systems with short range interactions, dynamical 
localization implies exponential decay of ground state correlations, up to an explicit correction. 
Second, the dynamical localization of random xy spin chains can be reduced to dynamical 
localization of an effective one-particle Hamiltonian. In particular, the isotropic xy chain in 
random exterior magnetic field is dynamically localized. 

1. Introduction 

One important reason for the popularity of quantum spin systems in theoretical and mathe- 
matical physics is that they can serve as relatively simple models to study interacting quantum 
many-body systems. In this setting, the physics of a single spin is governed by a finite- 
dimensional Hilbert space and hence (essentially) trivial. Therefore quantum spin systems 
are ideal models to focus directly and exclusively on physical phenomena arising from many- 
particle scattering. For the same reason one may want to study random quantum spin systems 
to better understand the effect of disorder on interacting quantum many-particle systems. This 
is our goal here. 

Naturally, there is a considerable amount of physics literature on such models and we will 
not attempt to survey these works here (but we give some references to papers closely related 
to our work in Section [4.11 below). Much less has been done on the rigorous mathematical 
level, where there was a short outburst of activity in the early 1990s, including the papers 
[30] , [H] , [31] and [3] . This came at the end of a decade which had seen rapid progress in the 
theory of the Anderson model and, in particular, the development of several methods which 
allowed to give rigorous localization proofs. The referenced works from the early 90s built on 
the results and strategies which arose from studying the Anderson model and found ways to 
apply some of these ideas to quantum spin systems with disorder. 

Our work here arose from a desire to return to these investigations from the point of 
view of progress made and new questions having been asked over the last two decades. One 
important development in this period was the increased attention to questions of quantum 
dynamics and, in particular, the development of new tools and refinement of existing tools to 
prove dynamical localization properties for the Anderson model and other types of random 
Schrodinger operators. A concrete goal of the present work is to discuss the meaning of 
dynamical localization in the context of quantum spin systems and to rigorously prove this 
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for at least one class of models. Let us say from the outset that our model of choice, the 
xy spin chain in a random exterior magnetic field, is particularly simple as its localization 
properties can be directly reduced to those of the Anderson model. A future goal should be 
to understand dynamical localization properties for other models of disordered quantum spin 
systems, including mult i- dimensional systems. 

Another reason for returning to the study of disordered quantum spin systems is re- 
cent progress in rigorously understanding localization for the Anderson model with electron- 
electron interactions. The papers [T5| [T6] and [5] showed how to extend the methods of 
multiscale analysis and fractional moments, respectively, to obtain localization results in suit- 
able regimes for a system of N interacting electrons in a random potential. These methods do 
not yet allow to study the thermodynamic limit of an electron gas in a random environment, 
i.e. an infinite volume limit in which the number of electrons is proportional to the volume. In 
the physics literature, e.g. [8], it has been suggested that localization will persist in the case 
of low electron density. For a chance to investigate this rigorously it is necessary to grasp the 
consequences of many-particle scattering effects in such systems. Characterizing localization 
in terms of the properties of the dynamics, rather than the spectrum, seems to be better suited 
to these many-body systems. Interacting Anderson models, due to the infinite-dimensional 
single particle dynamics, are much harder to analyze than random quantum spin systems. 
Nevertheless, by starting with the latter one can hope to get insights on how many-body 
localization can arise. 

The dynamics of short range spin systems satisfy a locality estimate commonly known as 
Lieb-Robinson bounds [33]. These estimates can be understood as upper bounds on the group 
velocity of spin waves or, alternatively, as bounds on information propagation in a spin system 
(a term inspired by quantum information theory, where quantum spin systems are thought 
of as models for interacting qubits). As reviewed in Section [2] below, Lieb-Robinson bounds 
can be proven for very general classes of quantum spin systems and have recently found much 
renewed interest and many applications. 

We will define a strong form of dynamical localization for a quantum spin system in terms 
of the validity of a zero-velocity Lieb-Robinson bound, e.g. Definition 12.21 Given this, our two 
main results can be summarized as follows: 

(i) If a zero-velocity Lieb-Robinson bound holds in a given quantum spin system, then 
- up to a logarithmic correction — the system satisfies exponential clustering, i.e. ground 

state correlations decay exponentially, independent of the system's size. This is the content 
of Theorem I2.3[ stated and proven in Section 12.31 below. Exponential clustering is known 
for quantum spin systems with a size-independent ground state energy gap. Thus a zero- 
velocity Lieb-Robinson bound can substitute for a ground state gap in the proof of exponential 
clustering. 

(ii) We will prove that the isotropic xy spin chain in an exterior random magnetic field 
satisfies a zero-velocity Lieb-Robinson bound after averaging over the disorder. This will be 
stated as Corollary 14.11 More generally, we review in Section 13.11 the well known reduction 
(e.g. [M]) of the anisotropic xy chain to an effective single particle Hamiltonian. We then show 
in Theorem 13.21 that dynamical localization of the single particle Hamiltonian implies a zero- 
velocity Lieb-Robinson bound for the anisotropic xy chain. In the isotropic case the single 
particle Hamiltonian is given by the Anderson model, where the corresponding dynamical 
localization property is well known. 
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In Section 14.21 we combine these two results to conclude exponential clustering for the 
isotropic random xy chain, see Theorem 14.21 This requires some modification of Theorem 12.31 
and its proof, due to the fact that Theorem 14.21 is not deterministic but refers to expected 
values of the relevant quantities. This result is closely related to earlier work in [30] where a 
direct proof of correlation decay for the isotropic random xy chain is given. In contrast to our 
result, the proof in [30] does not proceed via dynamical localization. We will further discuss 
this, as well as other related work in e.g. [29], within the concluding remarks of Section [5j 

Let us finish this introduction by stressing what we consider to be the most important novel 
contributions of our work: (i) We propose the new concept of a zero-velocity Lieb-Robinson 
bound as a means to describe dynamical localization of many-body quantum systems. Our 
definition, which we can verify holds for certain models, is a simplified version of the mobility 
gap proposed by Hastings in [21], see Section 12.31 for more discussion. In particular, we 
provide a general argument which shows that this form of dynamical localization implies 
exponential decay of ground state correlations, (ii) As a first example of a many-body system 
where a zero- velocity Lieb-Robinson bound can be rigorously established we study disordered 
xy chains. The core argument is to reduce the analysis, via a Jordan- Wigner transform, 
to dynamical localization for Anderson-type single particle Hamiltonians. One obstacle in 
proving our estimates, which we stress are uniform in time, is to overcome the non-locality of 
the Jordan- Wigner transform, and this is accomplished in Section 13.21 This argument goes 
beyond previously known non-rigorous results in the physics literature (see Section 14.11 for a 
discussion of the results in [TB"]). 

2. Lieb- Robinson Bounds and Correlation Decay 

The goal of this section is to prove an estimate on the decay of correlations in systems 
that satisfy a strong form of the Lieb-Robinson bound, see Definition 12.21 below. Although 
we state and prove this bound for deterministic systems, the most interesting models for 
which these bounds have been established are random, more on this in Sections |3] and HJ 
We begin this section by introducing a general class of quantum spin systems to which our 
result applies. Next, we briefly discuss Lieb-Robinson bounds and a result on exponential 
clustering for gapped systems. Finally, we define what we mean by a zero-velocity Lieb- 
Robinson bound and then state and prove our main result on correlation decay in systems 
which are dynamically localized in this sense, Theorem 12.31 below. 

2.1. Quantum Spin Systems. Our results on correlation decay apply to a very general class 
of mult i- dimensional quantum spin systems. In fact, while we choose the set of spin sites to 
be Z 17 , with some more notational effort we could equally well consider systems over much 
more general graphs. We feel, however, that this would not be beneficial to our presentation. 

A quantum spin system over r L v is defined as follows. To each m e Z", associate a finite 
dimensional Hilbert Space "H. m = C Um . l-i m is commonly referred to as the single-site Hilbert 
space, and the dimension n m > 2 is related to the spin J m at site m by n m = 2J m + 1, e.g., 
n m = 2 corresponds to spin J rn = 1/2. For any finite set iV C Z^, a composite Hilbert space 
and algebra of observables are defined by setting 

n N = (g) H m and A N = (g) B{H m ) , 

where B(H m ) denotes the bounded linear operators over Tt m , i.e., the set of n m x n m complex 
matrices. Due to the tensor product structure, it is clear that for any finite sets Nq C iV C Z^, 
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each observable A G An can be identified with an observable A' = A <g> ILyyvb e An- In this 
case, we regard An C An and thereby define inductively the algebra of all local observables 

*<4ioc = {J An 

where the union is taken over all finite subsets of 7L V . 

A model on such a quantum spin system is defined through an interaction. An interaction is 
a mapping $ : V (Z U ) -> A\ oc satisfying $(M)* = $(M) G A M for all finite M C Z". Corre- 
sponding to any interaction, there is an associated family of local Hamiltonians, parametrized 
by finite subsets N C Z", given by 

(2.1) H N =J2®(M). 

McN 

When the interaction is fixed, we will often drop the dependence of the local Hamiltonian on 
$. Since the sum in (12. ip is finite, Hn is self-adjoint. By the spectral theorem, there is a 
one-parameter group of automorphisms, which we denote by {rf*}t&., defined by setting 

if {A) = e itHN Ae- itHN for all A e A N . 

is called the Heisenberg dynamics or time evolution corresponding to Hn- 
The above general framework covers a wealth of heavily studied explicit models, such as 
the Heisenberg, Ising, XY and XXZ models, all describing spin 1/2 particles. An important 
mathematical model for higher spin is the spin 1 AKLT model pQ. For general background 
on the mathematical treatment of quantum spin systems as models of statistical mechanics 
we refer to [9] and, for a different perspective, [16]. 

As an example we provide the spin 1/2 Heisenberg model, probably the oldest and most 
prominent example. 

Example 2.1. A common model is the Heisenberg Hamiltonian. In this case, one takes as 
single-site Hilbert space 1-L m = C 2 for all m G 7L V . The interaction for this model is given by 

§< M \ = [ V ( a n a ™ + a n a m + K^m) if M = {n, m} and \n-m\ = 1, 
\ otherwise, 

where the real number \x is the parameter of the model, 

(2.2) I), 7), an, „■ = ( J \ ) 

are the Pauli-spin matrices, and for any w G {x, y, z}, a™ Am is the matrix a w <S> 1 with 
a w in the n-th factor. 

2.2. Some Prior Results. In 1972, Lieb and Robinson proved a locality estimate for a large 
class of quantum spin models [33]. Their observation can be described as follows. Consider 
a quantum spin system over and let J and K be finite, disjoint subsets of If . Take any 
finite subset N with JUK C N. It is clear from the tensor product structure of the observable 
algebras that for any A G Aj and B G Ar, [A ® ~&n\j, B ® ^n\k] = 0. To ease notation, we 
will often suppress these identities and simply regard A and B as observables in An whenever 
J,K<zN. For a large class of essentially short-range interactions, Lieb and Robinson proved 
a bound of the form 

(2.3) || [r^(A), B] (I < C(J,K)\\A\\\\B\\e-^ K) - v ^\ 
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where d(J, K) := min{|j — k\ : j G J, k G K}. 

This estimate shows that for times t with \t\ < d(J,K)/v the commutator remains expo- 
nentially small, and so disturbances do not propagate through the system, by the Heisenberg 
dynamics, arbitrarily fast. It is important to note that the numbers C(J, K), rj, and v are all 
independent of the volume N on which the dynamics is defined. The number v, which can 
be made explicit in terms of an appropriate norm on the interaction, is called a bound on the 
Lieb- Robinson velocity of the model under consideration. Recently there have been a number 
of generalizations of these Lieb- Robinson bounds J39J |26j [361 EL71 El EH ESI 0U H3] , and they 
have proven useful in a variety of applications [221 SOI [23l [2TI [38l [TTJ, [10l [7] . Many of the 
previously mentioned results were reviewed in e.g. [121 [251 HZ]- 

One of the first applications of the new Lieb- Robinson bounds was a proof of the Exponential 
Clustering theorem [321 [25]. Stated simply, exponential clustering describes the fact that 
ground-state expectations of gapped quantum spin systems decay exponentially in space. As 
one might imagine, the rate of exponential decay depends on the size of the gap, see [39] for an 
explicit estimate and [H] for an improvement. It is clear that local Hamiltonians corresponding 
to any interaction will be gapped. The important point in the proof of clustering is that if 
the model under consideration has a uniform positive lower bound on the gaps in the finite 
volume, then the rate of decay persists in the thermodynamic limit. 

Despite the fact that ground state expectations are time-independent, the proof of this 
clustering result uses Lieb- Robinson bounds in a crucial way. In fact, the proof requires a 
slightly stronger bound than the one claimed in (12. 3ft : 

|| [rf (A), 5] || < \t\C(J,K)\\A\\\\B\\e~ v{d{J > K) - vm , 

at least for small times, e.g., \t\ < 1. 

2.3. Correlation Decay in Dynamically Localized Systems. Our first result in this 
work, see Theorem 12 . 3 1 below, is a proof of a clustering-type bound under different assumptions 
than we described above. In words, we prove that a zero-velocity Lieb- Robins on bound implies 
exponential decay of correlations, up to a logarithmic correction. Interestingly, our method 
of proof, motivated by some recent observations by Hastings in [23], provides a decay rate 
independent of the size of the gap, and our bound depends on the gap only through the 
logarithmic correction. 

Before proceeding with this result, let us relate our work to the contents of [24]. There the 
concept of a mobility gap is proposed, which is considered a many-body version of dynamical 
localization (and inspired by the mobility edge concept for single particle systems). It is dis- 
cussed that systems with a mobility gap have properties similar to gapped systems, including 
decay of ground state correlations. They are also said to satisfy a higher-dimensional Lieb- 
Schultz-Mattis theorem and Hall conductance quantification. A zero-velocity Lieb-Robinson 
bound in our sense can be considered a special case, and particularly strong form, of a mobility 
gap. This allows for a simplified argument, although with methods similar to those of [21], in 
our derivation of correlation decay below. 

The goal of this section is to state and prove Theorem 12.31 below. We begin with the 
following definition, which provides a way to characterize dynamical localization of a quantum 
spin system. 

Definition 2.2. Let $ be an interaction on a quantum spin system over 7L V . We say that $ 
satisfies a zero-velocity Lieb-Robinson bound if there exists r] > with the following property: 
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Given any finite, disjoint subsets J,K C Z", there exists C(J,K) < oo such that 

(2.4) || [r t N (A), B] || < C(J, K) min[|t|, 1] \\A\\ \\B\\e** J & 

for all finite JVcZ" with J U K C iV, allt eM and all A £ Aj and B £ A K - 

Here it is important that the numbers rj and C(J,K) do not depend on the size of N on 
which the dynamics is defined. In applications one will usually be able to say more about the 
dependence of C(J,K) on J and K, see, e.g., Section [31 For bounds of the form (12. 4p to be 
useful, C(J, K) should only depend on simple geometric properties of J and K, such as their 
size or surface area. 

Also note here that IH/r^A), B]\\ = \\[A,T^ t (B)]\\, so that bounds such as (12.4p can be 
applied equally well to ||[A, r/ v (5)]||. 

Theorem 2.3. Let $ be an interaction on 7L V that satisfies a zero-velocity Lieb-Robinson 
bound. Let J and K be finite, disjoint subsets of 7L V and take a finite set NaT/ with 
J U K C N . Suppose the local Hamiltonian, H^, corresponding to $ satisfies miner (.f/jv) = 0. 
Let ipo = ipo(N) denote a normalized ground state, i.e., H^ipo = and \\ipo\\ — 1 ond let 
P is the orthogonal projection onto ker(ifjv). Then for all A G Aj and B G Ak with 
P Bip = P B*ij; = 0, the bound 



(2.5) |(^o,^o)| < 



1+^1*1 2-ln 



\A\\\\B\\e- ,ld{J > K) 



7T y ^TTf]d(J,K 

holds. Here t] and C(J,K) are as in Definition \2.2\ 7 = 7^ is the spectral gap, i.e. 

7 = mm(a(H N ) \ {0}) . 

Before we begin the proof of this theorem, some comments are in order. First, if the 
ground state of Hn is non-degenerate, then the additional assumption on the observable B is 
equivalent to Bip ) = 0. In this case, our result then demonstrates exponential decay of 
correlations \ {4>o, AB^ ) — (^0, Aip ){ip Q , Bifj Q ) \ in the distance of J and K (apply Theorem 12.31 
to B:= B-(^B^)1). 

Also, the assumption that has ground state energy is inessential and made exclusively 
for convenience in proofs. For more general Hm with ground state energy Eq one gets the 
same result by first applying Theorem 12.31 to '■= — Eq , which has the same ground 
state and Heisenberg evolution as H^. 

Next, for general systems and especially large sets N, we expect that 7 is quite small. In 
this case the meaning of ( 12. 5ft is 

l^o, ABtP )\ < ||A||||S|||ln7|e-("-) d W K ), 

in the usual sense that 77— can be 77 — e for and e > and < allows for an e-dependent (and 
(J, independent) constant factor. Clearly, the rate of decay is independent of the gap size; 
up to the logarithmic correction in 7. The usefulness of this estimate depends on the context 
of the application and, in particular, how rapidly 7 approaches zero for large volume. We will 
discuss this for the example of the isotropic random xy chain in Section |4] below. 

Finally, similarly to [39] and [26], see also [21], our method of proof employs a specific 
approximation technique. Fix finite sets K C TV C Z w . For any local observable B G Ak, 
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a > 0, and e > 0, a quasi-local approximation of B is defined by setting 



f 



-act? 



If K C N, then for each positive choice of a and e, the support of B(a, e) is N, however, the 
explicit gaussian kernel enables useful estimates. In fact, the following is a basic fact about 
Fourier transforms of gaussians which we will use in our proof. Let E e K, a > 0, and e > 0. 

1 f e iEt e- at2 1 f°° (»- b) 2 

(2.6) — / ^— dt = —= / e- ew e- L T^dw. 

27rz ./« t - 2 x /7ra / n 



This follows, e.g., from the identity 
(2.7) (e + zt)- 1 = / e~ w{e+lt) dw . 



oo 



Proof of Theorem \2.3[ We begin by observing that 

(2.8) (Vo, AB^o) = OAo, A (fl - S(a, e)) ^o) + (^o, B(a, e)A^) + <^ , [A, 5(«, e)]^ ) 

for any choices of positive a and e. An application of the spectral theorem and (12. 6p shows 
that for any x G Hn, 

r°° i r e itE e -at 2 

(x,B(a,e)ijj ) = — — ; — dtd(x,P E Bip ) 

Jo 2tt2 J r t-ie 



/ r-p= / e- ew e- J ^dwd{x : P E Bip ) 



where, for the last equality, we used the assumption that PoBip = 0. The integrand can be 
re-written as 



1 



°° _ (™-E) 2 

e ew e 



e dw H == / (e~ ew — 1) e ±^ dw 



=: l-fli(E,a)+i?2(£;,a,e). 
Since limsup e ^ Pl2(E, a, e) = and 2Ri(E, a) < e ^ for > 7, it is clear that 

/>oo 

lim sup I (t/j ,A(B -B(a,e))ip ) \ < / R\{E, a)d(A*tp , P E B^ ) 

e^O J~ 



< -e-£\\A*if> \\\\B1>o\\. 



Similarly, 

"oo 1 p „—itE „—at 2 



(Vo,5(a,e)AVo> = \ 7T-.\ — + : dtd(B*i> , P E Ai/; ) 

Jo 2m Jr t-ie 

R 3 (E,a,e)d(B*^ ,P E Aip ), 



where 

1 f°° _ {w+E) 2 1 _g2_ 

R 3 (E,a,e) = — -= I e ew e 4 <* dw < -e 4 <* . 

2Wna In 2 
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The bound 

1 7 2 

limsup|(^o,5(a,e)^ >| < o e ~^ ll^oll ll^oll , 

readily follows. 

To the final term, we apply the assumption of a zero- velocity Lieb-Robinson bound. Clearly, 
(2.9) \{^ ,[A,B(a,e)]^)\ < \\[A, B(a,e)]\\ <±-f Ml^Me'^ dt . 

We divide the integral above into three pieces. Using (12 .4ft . it is clear that 

2tt J\ t \<i \t\ TV 



and similarly, for any A > 1, 



2tt 

Lastly, 



- / iM.-'i < ^M)| M n| B || e -^) ln(A) . 

!vT A<|t|<A 1*1 ^ 



JL /■ II WW, a] 11^ < c«iO 

2jr./|,|>A |t| " 2tt AV a 

Now with the choice of 

(2 - 10) » = and 2A = ^- 

we have proven that 

|(^o,^o)| < \\A\\\\B\\e-^ (l + ^V^ (2 + ln^)) , 
as claimed. □ 



3. The xy -chain with disorder 

In this section we consider the anisotropic xy-chain with free boundary conditions, as in- 
troduced in jM]. We start by reviewing the diagonalization of this xy Hamiltonian with 
deterministic coefficients. In essence, this consists of using a Jordan- Wigner transform to re- 
duce the n-body xy Hamiltonian to a free Fermion system whose diagonalization is governed 
by the diagonalization of an effective one-body Hamiltonian. The latter takes the form of a 
2 x 2-block Jacobi matrix. 

Next, we consider the case where the coefficients in the xy chain are random variables. In 
this case the corresponding block- Jacobi matrix is a random operator closely related to the 
Anderson model. We demonstrate that if the block- Jacobi matrix is dynamically localized, 
then the Heisenberg dynamics corresponding to the xy Hamiltonian satisfies a zero-velocity 
Lieb-Robinson bound in average. 
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3.1. Diagonalizing the xy Chain. Consider three real- valued sequences {fij}, {7j}, and 
{z/j }. These parameters will represent the coupling strength, the anisotropy, and the external 
magnetic field respectively. We will assume fij 7^ for all j as otherwise the chain decom- 
poses into shorter pieces. To each integer n > 1, we will denote by -ff[i, n ] the finite volume, 
anisotropic xy Hamiltonian with free boundary conditions given by 

n— 1 n 

(3-1) #[i,n] = £>i[(l + 7iW + i + (1 - 7iW + i] + E • 

i=i j=i 

Here <7 X , a 2 ', and cr z are the Pauli matrices given in (12. 2p . and Hu n ] acts on the Hilbert space 
H[i, n ] — ®j=i ^- 2 - For notational convenience, we will write H n = H[i >n ] below. 

It is well-known, see e.g. [3l], that Hamiltonians of the form ( 13. ip can be diagonalized 
in terms of a system of operators satisfying canonical anti-commutation relations (CAR). To 
introduce useful notation, we briefly discuss this diagonalization procedure below. 

First, one introduces raising and lowering operators as 



a* = -(of + ia v A and a,- = -(of — ia y A for each 1 < j < n. 
j o ■> ■> 2 — — 



1 , „ . „s . 1 

Using the relations 

n x x _i_ y y _ j 2 (a*a j+ i + a* +1 a,,) if +, , z _ „ * fl 

± <W+i " j 2 (a,a i+1 + a* +l a*) if -, aM ^ " 2a ^ 1 ' 

one sees that 

n— 1 n 

(3.2) = 2 ^ fij[a*a j+1 + a* +1 aj + 7i(ojOj+i + + E ^( 2a j a i ~ J )- 

i=i j=i 

Since the a-operators preserve the local structure, they do not satisfy anti-commutation rela- 
tions. To break this locality, one introduces new variables via the Jordan- Wigner transforma- 
tion. Set 

(3.3) ci = a\ and Cj = o\- ■ • o z ^_ x aj for all 2 < j < n. 
A short calculation shows that 

c^-Cj = Qjjdj) CjCjj r \ = djCLjj^i^ and CjCj_^i = dj(Xjj r \ , 
in which case we find that 

n— 1 n 

(3.4) H n = 2j2 V-i \~ c *j c i+i - c i+i c i + Tj ( c i c i+i + c i+i c j )] + E "i ( 2c i c i ~ n ) ■ 

In contrast to the a-operators, the c-operators do satisfy the CAR. In fact, one easily 
calculates that for all 1 < j, k < n, 

[ ' { Cj , Cfc } = {c*,c*}=0, 

which are the canonical ant i- commutation relations. It is convenient, and equivalent, to 
express these CAR through a vector-valued formalism. Let C = (ci, . . . , c n , c*, . . . , c*)*, then 
it is clear that the c-operators satisfy (13 ,5p if and only if 

(3.6) CC* + J(CC*fj = 1, 
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where 



J 



II 
11 



J 1 - 



Here and in the following A 1 is the (non-hermitean) transpose of a (not necessarily square) 
matrix A, i.e. its matrix-elements are (A^jk = Akj- We will use this notation for matrices A 
with scalar entries as well as for the case of operator-valued entries. Technically, C, J as well 
as many of the quantities introduced below depend on n, just as the Hamiltonian H n , but we 
will drop this from our notation to avoid over-subscripting. 

We find it useful to express the xy Hamiltonian H n in terms of C. Symmetrizing (I3.4p using 
(I3.5p . it is clear that 

n—l n 
3=1 3=1 



n-l 



J2^3 [- 



■CjCj+i + Cj+iCj - C j+1 Cj + CjC 



■j+ 



3=1 



n-l 



* C* 



(3.7) 



3=1 

C*MC, 



J + 

n 

C J C *3+l] +Y U j( C *3 C 3 ~ C 3 C *j) 
3=1 



where M is the 2 x 2-block matrix 

(3.8) M = 

with A and B Jacobi matrices satisfying 

/ v\ —Mi 



A B 
-B -A 



(3.9) 



and 



(3.10) 



.4 



\ 



-fJ>i 



V 

/ -/ii7i 
Mi7i 



—fJ>n-l 
—fJ"n-l v n ) 



B 



\ 



-f^n-lln-1 

Observe that A* = A 1 = A, B* = B l = -B, and thus M* = M* = M. 

We complete the diagonalization of H n by applying a Bogoliubov-type transformation. To 
see this, set S = A + B and denote by < Ai < A2 < . . . < A n the singular values of S 
counted with multiplicity. The singular value decomposition of S gives real-valued orthogonal 
matrices U and V such that 



(3.11) 



USV = U{A + B)V l = A, 
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where A = diag(Ai, . . . , A n ). In this case, 

A = A* = v&U* = V(A - B)U\ 

Set 

w= l -( v + u v ~ u 

2\V-U V + U 

and note that W is an orthogonal matrix. This can be checked by directly verifying that 
WW 1 = t, or, alternatively, by noting that 

an orthogonal matrix. A short calculation shows that W diagonalizes M, i.e., 
(3.12) W ( A u B A )w f 



B -A I \ -A 

Finally, we define 6-operators by setting 

(3.13) B = WC. 

It is easy to check that B has the form 

# = (pi, . . . ,b n ,b\, . . . , 6*)*. 

Moreover, since \J , W] = 0, it is clear that 

BB* + J{BB*YJ = W (CC* + J(CC*YJ) W l = 1 

i.e., the fe-operators satisfy CAR as well. 
From (13.71) . we conclude that 

H„ = CMC = B* ( ^ ° i )B 



-A 



n 



(3.14) = 2j2^b*b j -E^t, 

3=1 

where = E"=i ^i- 

This means that when written in terms of the 6-operators the xy Hamiltonian takes the form 
of a free Fermion system. This allows to explicitly diagonalize H n (in as far as the singular 
value decomposition (13. lip is explicit), which we will exploit further in the next section for 
the case of the isotropic xy chain. 

We close this section by observing that (I3.14p makes it easy to compute the time evolution 
of the 6-operators. In fact, for any 1 < j < n, 

r™ (b J )= e~ 2itx i bj and rf (b*) = e 2UX > b* , 

where r t n is the Heisenberg dynamics corresponding to H n . This follows, for example, by 
observing that, due to (13.141) and CAR, [H n ,bj] = —2Xjbj and thus 

jT^bj) = ir?([H n , bj]) = -21^%) and r^) = b 3 . 
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In vector form, this can be written as 

/ -2UA q \ 

<{B) =[ e Q J itA J B. 

More is true. Since W is a matrix of scalars, it is clear that 

7f(C) = W 1 t?{B) = e~ 2itM C, 

or, componentwise, 

n n 

(3.15) 7f(c,) = £ M hk (1t)c k + Mj, n+k (2t)cl 



k=l k=l 



where 



M j>k (t) = (e- iMt ). k . 
The above facts will be useful in the following sections. 

Remark: The specific structure of M suggests the change of basis (ei, t2, . . . , &2n) to 
(e±, e n+ i, 62, e n+ 2, . . . ,e n , t2n) in representing M. This leads to 

/ viJ -/ii5'(7i) \ 



(3.16) M = M 



Here 



-^n-lS^n-l] 
-/i„_lS'(7 n _i)* V n J 



1 \ or \ _ ( 1 7 

-7 -1 



J = o -i • 



In this view M takes the form of a generalized tight-binding Hamiltonian with a (sign- 
indefinite) potential term generated by the magnetic field (uj) and non-standard hopping 
terms created by the 2 x 2-matrix-valued off-diagonal terms — /ijS , (7 J ) and their adjoints. 
Mathematically, (I3.16P provides the possibility to investigate spectral properties of M with a 
transfer matrix formalism, although of higher order than in the case of standard tri-diagonal 
Jacobi matrices. 



3.2. From Dynamical Localization to Zero Velocity Bounds. In the previous sub- 
section, we demonstrated how to diagonalize a general xy Hamiltonian. The goal of this 
subsection is to prove that dynamical localization of the random matrix M = associated 
with H n via (13. 7p implies a zero- velocity Lieb-Robinson bound for H n . 

We begin by discussing what it means for to be dynamically localized. Let the infinite 
sequences {fij}, {7j}, and {vj}, which define an xy model and thus the matrices M^ n \ see 
Section [37T] correspond to sequences of random variables over a probability space 
some concrete examples are provided in the following section. Thus the are self-adjoint 
random matrices and we will define dynamical localization in terms of bounds on the expected 
values of the matrix elements of their time evolution. In the following E(X) = j n X dP denotes 
the expectation of a random variable X on Q with respect to P. 
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Definition 3.1. We say the matrices M^ 1 ' are dynamically localized if there exist numbers 
C > and rj > such that for any integers j, k, n > 1 with j, k G [1, n], 

(3.17) E (sup \M$(t)\ + sup \M^ +k (t)\) < Ce~^ . 

Here M^ n k \t) = ( e ~ lM{n)t \ ; where is the matrix associated to H n via ^3. 7| ). 

Dynamical localization can be equivalently expressed in terms of the rearranged matrix 
M = AfW from fl3TT6|) . 



?j|fc-j| 



E I sup 

\tm 

for all n and j, fc G [1, 2n] with slightly modified constants C and rf. 

It is important to note that the numbers C > and r/ > in Definition 13.11 are independent 
of n. In this case, the bound in (13.171) is uniform with respect to arbitrary finite volumes. The 
next section describes a class of xy models which can be shown to satisfy Definition 13.11 

The main result of this subsection is that dynamical localization of implies a zero- 
velocity Lieb-Robinson bound after averages over the disorder are taken. 

Theorem 3.2. Assume that is dynamically localized in the sense of Definition \3.1\ 

There are numbers C > and 77 > such that for any integers 1 < j < k and any n > k, the 
bound 

(3.18) e(su V \\[t?(A),B}\\) <C"||A||||£||e- 

\tm J 

holds for all A G Aj and B G A[k, n ] ■ 

In fact, with C and rj from ( 3.11 ), we may choose the same exponent rj and 

Proof. Fix n > k. Let us first examine the case that A = Cj. Using (I3.15p . it is clear that 

n n 

[7f (c,-), B] = M hk >(2t) [c k ,, B] + J2 M hn+k> (2t) [<£/, B] , 

k'=k k'=k 

where we have used that B G A[k, n ], cy, c* k , G A[i : k'\ an d thus [ck>, B] = [c* k ,, B] — for k' < k. 
The bound (I3.17P implies that 

(3.20) E (sup\\[rl l ( Cj ),B]\\ \ < 4C\\B\\ ^ e^'^ < 

^ * ' k'=k l-e " 

where we have used that \\[A, B}\\ < 2||A||||£>|| holds in general. By taking adjoints, the case 
of A = c* follows with an identical estimate. 

Now suppose that A = aj. By (13. 3p we have aj = erf . . . o~ z -_ x Cy Using the automorphism 
property of r t n and the Leibnitz rule, we have that 

[TfCay), B] = Tf K) • ■•7?(a' j _ 1 )[7?{c j ), B] + [Tf (a?) • •■r^_ 1 ), B]t?{c,). 
Clearly then 

(3.21) E (sup || [r^ajlB] ||) < ^M- e -v(k-j) + E ( sup || . . . B ] 
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where we have used (13.201) . We now further expand the second term above. 
For any 1 < £ < k, it is convenient to define 



C(£,5)=Efsup||[r t >f)-..r t >l),5]|| 
We claim that 

(3.22) C(£, B) < C(£ — 1,B) + ^£Ml e -v(k-e) _ 

To see this, note that another application of Leibnitz shows 

(3.23) \\[T?{al)---T?{al),B]\\ < ||[t? (of), B}\\ + || [ T » (of) • • •t'K-i), B] || . 
Moreover, for any 

erf = 2a^a£ — 11 = 2c|q — 11 , 

and so 

(3.24) [Tfto), £] = 2 [t?(c2), B] r«( Q ) + 2r»(cJ) [r^( Q ), B] . 

It is now clear that fl3T22|) follows from fl3~23|) . fl3~24j) . and fl3~20|) . Additionally, the quantity 
corresponding to I = 1 satisfies 

c( l ; m<f^M e -^-i) ) 
v 7 _ 1 - e-n 

as is clear from (I3.24p . 

As a result, our estimate now follows by iteration. From (13.2 ip . we have that 

E(sup\\[T?( aj ),B}\\\ < l3He-^ + C(j-l,B) 



< fgM e -^(H) f 1 + e ~ vu ~ e) 

4?C\\B\\ ^ 
~ (1-e-") 2 

Again, the case of A = a* follows, with the same estimate as above, by taking adjoints. 
Using Leibnitz, the above bound yields estimates for both A = a*aj and A = Oja*, increasing 
the bound only by 2. Since the collection {aj *} constitutes the canonical basis 

for Aj, we get ( 13.1 8ft with the explicit constants ( I3.19p . using the simple fact that the norm 
of a matrix bounds the absolute values of all of its entries. A slightly more refined argument 
would allow to lower the constant 96 in (I3.19P a bit. □ 

A number of comments about Theorem 13.21 are in order. 

Remark 3.3. Theorem \3. 6 A applies to more general observables. For example, let J C Z be 
a finite set. To each j G J and any 1 < n,m < 2, denote by ej(n,m) G Aj a matrix unit, 
i.e., a matrix with 1 in the (n,m) entry and all other entries equal to 0. Any A G Aj can be 
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expanded in terms of products of matrix units; there will be 2 2 ' J ' such terms. Here \J\ denotes 
the cardinality of J. Using Leibnitz to expand the product, it is clear that 



holds for all A £ Aj and B £ A\k, n ] if max[j : j £ J] < k < n. Here d(k, J) = min[|&; — j 



The bound ( I3.25P is useful in applications where J is kept fixed as n varies. However, when 
applying ( 13.251) to cases such as A £ A[ij], B £ A[k, n ], j < k, the constant on the right grows 
exponentially in j, while one would hope to have bounds exponentially small in \k — j\ and 
constants uniform in n and 1 < j < k < n. 

The next result, which we state as a corollary, achieves a bound of this type, holding for 
a larger class of observables than Theorem 13.21 Moreover, this result provides better bounds 
for small times, e.g. \t\ < 1. The drawback of this result, and the method to prove it, is that 
one pays a penalty for large times. 

Corollary 3.4. Consider an xy model that is dynamically localized in the sense of Defi- 
nition \3.1[ If all of the random variables are uniformly bounded, then there exist num- 
bers c > and 77 > such that for any finite set J G N and any integers n, k satisfying 
max{j : j £ J} < k < n the estimate 



Here / depends on n, A, and B, but we will suppress this in our notation. Calculate: 



Above we have denoted by Hj those terms in H n that have non-trivial commutator with 
A, and the final equality follows from the Jacobi identity. This proves that / satisfies the 
differential equation: 



(3.25) 





r]\k-j\ 



fit) 



i[r?([H n ,A]),B] 
i[T?{[Hj,A]),B] 
i][<r?(Hj),<T?(A)],B] 

-i[[T?(A),B],T?(Hj)]-i[[B,T?(Hj)],T?(A)\ . 



f'(t) = -i [f(t), t?(Hj)} - i [[B, t?(Hj)}, t?{A)] . 



Since the first term is norm-preserving and /(0) = 0, we get 



(3.26) M{A\B\\\ = \\f(t)\\<2\\A\\ f \\K(Hj) : B]\\ ds , 




Hj = 5>[(1 + 7iW+i + (1 " 7i)*M+i] + £ v i°i ' 
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where J = [ J U ( J — 1)] H [l,n]. For any w G {x,?/}, the estimate 

E (II fc^X+i), B]\\)< <7'(1 + e")||5|| e -"M , 

follows from (13. 18[) by Leibnitz; it holds for each j G J and is uniform in s > 0. Using fj3.26j) 
then, we find that 

1*1 



E(||[7?(A),5]||) < 2\\A\\ E (|| [t™(Hj), B}\\) ds 



o 

< AC\\A\\^C'{\ + e 1, )\\B\\e- r '^- i) 



t\ + 



+2C\\A\\J2 c '\\ B \\ e ~ ri{k ~ j) 
\ B \\ e -vd(k,J) 



t\ 

< c\t\\\ ;;;; "' '\ 

where C = max[|^|(l + \^\)\, C = max[\uj\] and c := (4(7<7'(1 + e") + 2(7(70/(1 - e""). □ 

4. The isotropic random xy chain 

In this section we will apply the above results to the isotropic xy chain in an exterior random 
magnetic field. This corresponds to the special case jj = and \ij — \i ^ for all j in (13. ip . 
giving the Hamiltonian 



(4.1) 



We will assume that 



i=i j=i 



(A) The coupling parameters uj, j G N, are i.i.d. random variables, whose common distribu- 
tion Po is absolutely continuous with bounded and compactly supported density p. 

In this case one may choose (fi, JA,F) as the infinite product space ,- eN (K, B, P ) with 
the Borel algebra B. 

4.1. Dynamical Localization. For the isotropic xy chain the matrix M in (I3.8P becomes 
block-diagonal, 

' A 
-A 



(4.2) 
where 



M 



( vi -(i 



A 



\ 



\ -fj, v n ) 

is the restriction of the classical Anderson model to [l,n]. Dynamical localization ( 13. 17ft for 
M becomes equivalent to dynamical localization for the Anderson model in the form 



(4.3) 



E(sup|(e- lA ') Jifc 
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for constants C > and r/ > uniformly in n G N and 1 < j, k < n. Under the above 
assumption on the distribution of this is well known. It has been proven by two different 
methods: 

(a) The Kunz-Souillard method goes back to [32], see [18] and [32] for detailed presentations. 
We remark here that these sources provide a proof of (14. 3 j) for the Anderson model over [— n, n] 
(instead of [l,n]) and for (j, k) = (0,k), < k < n. The method can easily be adjusted to 
handle general pairs (j, k). 

(b) The fractional moments method was initially developed by Aizenman and Molchanov in 
[I] and is first shown to yield dynamical localization in [2] . Detailed arguments leading to (14. 3[) 
can be found by combining the exposition in Section 6 of [IS] with fractional moment bounds 
for the one-dimensional Green function as provided by Proposition A.l in [35] or Theorem 4.1 
in [21]. 

In summary, we have the following corollary of Theorem 13.21 

Corollary 4.1. Suppose that the magnetic field in the isotropic xy chain OTTp satisfies (A). 
Then the Heisenberg evolution 

t 71 (A^ — e^^ n ' lso ^ Ae~ ^^ n ' is °^ 
satisfies A3.18\) . a zero velocity Lieb-Robinson bound in disorder average. 

Remark: Corollary 14.11 improves on a result stated in [T3] , where it was argued that the 
isotropic xy-chain in random magnetic field satisfies the bound 

(4.4) \\[rl l (A),B]\\ < cn 2 \t\e- vlj - kl 

for some 77 > and A G Aj, B G Ak- This was interpreted as a logarithmic light cone for 
information propagation in the spin chain. As in our approach, [T3] uses the Jordan- Wigner 
transform to reduce this to localization properties for the Anderson model. In the present 
work we use suitable rigorous versions of the latter to show that, after disorder averaging on 
the left hand side of (I4.4p . one gets strict dynamical localization in the sense that the right 
hand side neither grows in time nor in the size n of the chain. 

We will not attempt here to survey other results on disordered spin systems in the physics 
literature, but we mention several more works on closely related models: Numerical results 
on the dynamics as well as on spin-spin correlations of the Heisenberg XXZ chain in random 
magnetic field are presented in [35] . In [T2] quantum spin chains with temporally fluctuating 
exterior fields are studied and it is discussed that these lead to dynamical phenomena different 
from those observed for models with static disorder. Finally, a very comprehensive discussion 
of the physics of different types of spin 1/2 chains (in particular XY and XXZ) with random 
exchange couplings (e.g. the parameters fij and rjj in the case of the XY model) can be found 
in [20]. 

4.2. Decay of ground state correlations. Corollary 14.11 allows, at least in principle, to 
apply Theorem l2.3l to the isotropic xy chain in random magnetic field and conclude exponential 
decay of ground state correlations. This is not completely straightforward as there are two 
additional issues which need to be addressed: 

(i) Theorem 12.31 is deterministic in nature, saying that a zero- velocity Lieb-Robinson bound 
for a fixed quantum spin Hamiltonian implies exponential decay of ground state correlations for 
this Hamiltonian. Corollary 14. 1[ however, provides the zero- velocity LR bound in expectation. 
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To get a probabilistic version of correlation decay from this we will have to slightly adapt the 
proof of Theorem 12.31 

(ii) We need information on the gap size 7 for the random xy chain to control the logarithmic 
correction in (12. 51) . In doing so we will have to take into account that 7 itself is a random 
variable. To achieve the second goal, we start by providing some more background on the 
explicit diagonalization of the free Fermion system (I3.14p . 

As we are interested in the isotropic xy chain and thus M as in (14. 2 j) , we can simply choose 



W 



U 

u 



in (I3.12p where U is orthogonal and diagonalizes the Anderson model, 

tfAtf* = A = diag(Ai,...,A n ). 

Note that this choice is slightly different from what was done in Section 13.11 as Ai < . . . < 
A n are now the eigenvalues of the Anderson model, which are not necessarily non-negative. 
Proceeding very similar to above, in fact a bit simpler, we define 















V K ) 




V c n / 



see that bj, j = 1, . . . , n, satisfy the CAR and allow to represent the isotropic xy chain as the 
free Fermion system 

n 

H n ,- lso = 2j2W*b j -E( n H, 
i=i 

where, as before, = Y^j=i -\r The bj and bj are interpreted as creation and annihilation 
operators. The following facts are well known, see e.g. [4"tj] . 

• The operators b*bj, j = l,...,n, are pairwise commuting orthogonal projections, 
meaning that they can be simultaneously diagonalized. 

• The intersection of the kernels of b*bj is one-dimensional. Thus it is spanned by 
an essentially unique normalized vector Q (the vacuum vector), so that b*bjQ = 0, 
j = 1, . . . ,n. 

• For each a = («i, . . . , a n ) G {0, 1}™ use successive creation operators to define 

^ a = (b*) a tt := ■ ■ ■ iPl)^ 

which are an orthonormal basis of Hu^n] consisting of common eigenvectors for the 



0, if atj = 0, 



Clearly, the ip a are also eigenvectors for H, 



n,iso) 




H n , iso ip a = ( 2 x j - E(n) |i=| EM ^ 
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This shows that the ground state energy of H n is E = — Y^=i I A/ 1 anc ^ that -^o * s non ~ 
degenerate if and only if Xj ^ for all j. In this case the ground state is ip a (o), where 



(4.6) a 



(o) _ J 1, if Xj < 0, 
0, if Xj > 0, 



3 ' 

and the gap 7 between the ground state energy E and the first excited energy Ei is 

(4.7) 7 = J E 1 - J B = 2min|A j | = 2 dist(0, cr(A)). 

j 

Thus probabilistic lower bounds for the gap size derive from the following Wegner estimate 
for the Anderson model, which holds under our assumptions, see e.g. Theorem 5.23 in [2H] for 
a detailed proof. It guarantees the existence of a constant CV < 00 such that 

(4.8) P(dist(£, a {A)) < e) < C w \e\n 
uniformly in n G N, e > and £gR. For E = this implies by (14. 7p that 

(4.9) P(7 > 2e) > 1 -C w |e|n. 

From (14.81) we also see that P(0 G cr(A)) = and thus that the ground state of H n ^ so is almost 
surely non-degenerate. By ipo we denote this almost sure essentially unique and normalized 
ground state (the choice of phase is irrelevant for the following result). 

Theorem 4.2. Under assumption (A) there exist C < 00 and rj > such that 

(4.10) E(KVo,AB^ > - (^0,^0) (^0,^0) I) < C||A||||5||ne-^-*l 
/or a// A G A/, .B G ^4[fe,n] and all 1 < j < k < n. 

Proof. We closely mimic the proof of Theorem 12.31 with the main difference being that the 
zero velocity Lieb- Robinson bound holds in expectation (i.e. in the form (13.181) ). and not 
deterministically as assumed in Theorem 12. 3[ 

As the ground state of H n .\ so is almost surely non-degenerate we may assume by the remarks 
after Theorem 12.31 that (ipo, Bipo) = and, after an energy shift, that Eq = 0. We thus need 

to bound E(|(V>o,ABV>o)|)- 

Using Corollary 14.11 combined with an application of Corollary 13.41 to H niso yields the 
existence of C > and rj > such that 

(4.11) E(||[rr(A),S]||)<C'min[|t| ) l]||A||p|| e -^l 

for all A G Aj, B G A\k,n] an d 1 < j < k < n. 

Decompose Q into good and bad events according to 

n g ■= {v = (uj) g n : 7 > e-" n }, n b ■= n \ n g . 

By fT4T9l we have 

(4.12) P(0 6 ) < -C^ne""". 
Accordingly we decompose 

(4.13) E (\(ip ,AB*j; ) |) = E( X n g ■ \(^o,AB^ )\) +E( X n 6 ■ |^o,AB^ )|), 
where and xn,, denote the characteristic functions of Q g and fi&. 
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With ( I4.12p we can bound the second term by 
(4.14) 

E( X n b ■ \(ih,ABM) < \\A\\\\B\\¥(n b ) < ±C W \\A\\ \\B\\ n e^ n < ±C W \\A\\ \\B\\ n 

The first term on the right hand side of (14. 13j) is dealt with in complete analogy with the 
proof of Theorem 12 .3[ starting with the three-way split in (12. 8p . Due to the restriction to Q g we 
have 7 > e~ vn . Thus the proper choice of a corresponding to (12.1 Q[) is a := e~ 2vn / (4rj\k 
We again set A := ^Jtt /Aa. 

Proceeding as in the proof of Theorem 12.31 we get 

(4-15) E (xn g | {ipo, A(B - B(a, e))^) I) < ^\\A\\\\B\\e^\ 

and 

(4.16) E( X aMo, B(a,e)Atff )\) < ^||A|||| J B|| e -«. 
On the remaining term we use Fubini to bound, according to (12.91) . 

E(xn 9 |(Vo,[A^(«,£)^o)|) < E(\\[A,B(a,e))\\ 

As above, the t-integration is split into the the three regions |£|<1,1<|£|<A and \t\ > A. 
Proceeding as above, using (I4.1ip . one gets three bounds which combine to 

(4.17) E( X n a \(i(>oM,B(<x,e)]M) < + lnA) e -«. 

We finally observe that In A = r\n + ^ ln(7r?7| A; — The bounds (I4.15p . (I4.16P and (I4.17P 
combine into 

E ( X n g ■ I (^o, ABii> )\) < C\\A\\ \\B\\(1 + n + In \k - j\) e -^. 
Together with (I4.13P and (14.141) this gives (I4.10p . where rj' can be any number less than 7]. 

□ 

Remark: We note that Theorem 14.21 is trivial for the case of strong magnetic field in the 
sense that either suppP C [2\fi\, oo) or suppP C (—00, — 2|/i|]. In this case the ground state 
-00 of H nt [ SO has all spins down (or all spins up, respectively). Thus it is a product state with 
trivial correlations. 

To see this, say for supp Po C [2|/i|, 00), note that in this case all eigenvalues of the Anderson 
model A are positive, < Ai < . . . < A n . Thus, by (14.61) . the non-degenerate ground state of 
Hr, i SO coincides with the vacuum Q and 



i=i 3=1 



But one also verifies directly that the all spins down product state is an eigenvector of H n ^ so 
to eigenvalue — ■ Uj. 

Theorem 14.21 is non-trivial for the case where lies in the bulk of the spectrum of the 
Anderson model and thus the ground state of H n ^ so becomes entangled. 
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5. Concluding Remarks 

5.1. The work by Klein and Perez [30]. A result on exponential decay of ground state 
correlations for the isotropic x?/-chain ( 14.1 j) closely related to Theorem 14.21 was previously 
proven by Klein and Perez in [3D]. They consider the special case of correlations between the 
raising and lowering operators a* and and show that there exists an r/ > such that 

(5.1) sup 1(^0,^^0)1 <C v e-^- k \ 

n 

for all j and k and almost every choice of the magnetic field v = (vj)- Their proof is also 
based on the Jordan- Wigner transform and then uses Wick's Theorem to expand the ground 
state correlations (-00, aj, a^a) in terms of the Fermi two-point function (ip Q , c*Ckipo)- The 
latter is given by matrix elements of a spectral projection of the underlying Anderson model. 
Exponential decay is concluded using localization results for the one-dimensional Anderson 
model found via the multiscale analysis approach. Compared to our approach this has the 
advantage of allowing for general distributions of the i.i.d. random variables (vj), including 
the most singular case of Bernoulli variables. 

However, the results of [30] and our results on the xy chain are not directly compatible in 
that Klein and Perez get correlation decay with probability one, while we get our decay bound 
after averaging over the disorder. In principle, expectation bounds such as (I4.10p imply similar 
almost sure bounds with a reduced exponent by a standard summation argument. However, 
this would lead to a slight volume dependence in the almost sure bounds, due to the factor 
n on the right hand side of (I4.10p but also the fact that we work in finite volume and thus 
can't use translation invariance of the random variables considered here. Thus the almost 
sure bounds obtained from (I4.10p . while holding for larger classes of observables A, B, would 
be slightly weaker than the ones found in [30J. On the other hand, almost sure bounds such 
as (15.11) do not imply bounds on averages, due to the ^-dependence of the constant in (15.11) . 

The most important new aspect of our approach, as compared to the direct proof of cor- 
relation decay in [3D], is that we prove dynamical localization in the form of Corollary 14.11 
and then deduce correlation decay as a model independent consequence of this. This latter 
fact can be seen as somewhat reminiscent of a well-known result from the theory of random 
Schrodinger operators, saying that suitable forms of dynamical localization imply spectral 
localization (e.g. pure point spectrum with exponentially decaying eigenfunctions) . 

5.2. More general random block operators. Our strongest results on localization proper- 
ties for random spin systems, e.g. Corollary 14. II and Theorem 14. 2\ are restricted to the isotropic 
xy chain (14.11) . as in this case we can readily refer to the underlying dynamical localization 
bound (14.31) for the Anderson model. Extending these results to the anisotropic xy chain (13. ip 
would require to prove dynamical localization (I3.17P for more general random block operators 
M of the type (13.81) . (13.91) . (13.101) . under suitable assumptions on the random parameters Vj 
and Hj. Even for the case of constant [ij and i.i.d. Vj this does not seem to be known. A diffi- 
culty which prevents a straightforward extension of existing methods is the non-monotonicity 
of M in the random paramaters (the quadratic form of M is not monotone in the Uj, as each 
appears twice, once with positive and once with negative sign). Also interesting is the case 
of random interaction strength fij of neighboring spins (and, for example, vanishing or con- 
stant magnetic field), leading to an Anderson-type operator M with off-diagonal randomness. 
While the latter type of randomness has been studied in the literature, we are not aware of a 
result like (13.171) in this case. 
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A further generalization of the anisotropic xy chain ( 13. ip is found by allowing interactions 
of more than just next neighbors. By the same argument as in Section I3TT1 this can be reduced 
to diagonalizing a block operator M as in ( 13. 8ft with more general band matrices A and B. 

We also note that random operators of the slightly different form 



where both, A and £>, are self-adjoint have been studied in [29]. Here the authors were 

motivated by physical applications in the BCS-theory of superconductivity. For these models 
of random block operators as well, questions relating to localization remain widely open and 
will require additional work. 
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